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Stochastic Field evolution of disoriented chiral condensates 
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I present a summary of recent work Jl[ where we describe the time-evolution of a region of disoriented chiral 
condensate via Langevin field equations for the linear a model. We analyze the model in equilibrium, paying 
attention to subtracting ultraviolet divergent classical terms and replacing them by their finite quantum counter- 
parts. We use results from lattice gauge theory and chiral perturbation theory to fix nonuniversal constants. The 
result is a ultraviolet cutoff independent theory that reproduces quantitatively the expected equilibrium behavior 
of pion and a quantum fields. We also estimate the viscosity r}(T), which controls the dynamical timescale in the 
Langevin equation, so that the near equilibrium dynamical response agrees with theoretical expectations. 
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The relativistic heavy ion collider (RHIC) is 
currently smashing nuclei together to produce 
the most extreme states of nuclear matter ever 
achieved in a controlled environment. Their cool- 
ing dynamics may then reveal collective proper- 
ties of hot nuclear matter viz. its equation of 
state and the nature of the confining and chiral 
symmetry breaking transitions. 

Extracting this information requires dynami- 
cal extrapolation from the final products that hit 
particle detectors, across possible critical points, 
up to the initial hot phase. This is practically 
impossible in general but the theory of critical 
dynamics, a byproduct of the renormalization 
group, tells us how it may be achievable by mod- 
eling the dynamics of light degrees of freedom as 
stochastic (classical) fields. In nuclear Physics 
these are the pion and a fields and the associated 
equations of motion are stochastic versions of the 
Euler-Lagrange equations for the linear u-model. 

The simplest Langevin prescription is to add 
stochastic sources £ a (x, i) and dissipation terms 
to the classical field equations for <fi a = (a, 7?) 



dt 2 



V 2 a + A i 
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*I thank Krishna Raj agonal and Jim Steele for a stimulat- 
ing collaboration. See [M for a full account of the results 
summarized here. 



The stochastic fields £ a are taken to obey the 
fluctuation-dissipation relations 



0/)> = 2Vab(x)TS 4 (x-y) 



(2) 



with (£ a (x)) = 0. Here, we shall always choose 
Vab(x) = rj(T)8 a b- This dynamics reduces to the 
classical microcanonical relativistic evolution for 
the pions and sigma as r\ — > 0. 

In the limit of vanishing masses the long wave- 
length modes of the fields are effectively over- 
damped, and suffer critical slowing down [Q. 
Then our equations correspond to Model A, in 
the classification of dynamical stochastic theories 
of Hohenberg and Halperin || . Although Model 
A will be adequate for our purposes, a complete 
analysis of the universal 0(4) dynamics in the 
chiral limit of QCD requires coupling the order 
parameter (which is not conserved) to other con- 
served quantities in the theory, as in Hohenberg 
and Halperin's Model G §. 

The long time Langevin dynamics takes the 
system to its (classical) canonical ensemble, at 
temperature T. This results in divergences of cer- 
tain field expectation values as the lattice spacing 
a s is taken to zero. To make physical sense of the 
model in this limit requires a T dependent renor- 
malization scheme. This is achieved by making 
the replacement in the bare Eq. (1), 
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This removes linear and logarithmic divergences 
in the self energy of pion and a fields. 
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Figure 1. The order parameter {<fi} vs. tempera- 
ture T, for H = (circles) and H nonzero (trian- 
gles). Open and filled symbols show results ob- 
tained with a s = 1 on a 26 3 lattice and a s = 0.4 
on a 64 3 lattice respectively, with counterterms 
as described in the text. The line shows the pre- 
diction from chiral perturbation theory, Eq. (^). 

In Fig. [l], we required that the order parame- 
ter (</>} agrees with the expectations from chiral 
perturbation theory at low temperature M, 



{ct>)=U[l- 



T 2 
12/ 2 



(4) 



In particular, we have enforced the absence of 
linear T-dependence at small T by the addition 
of a further mass counterterm. This countert- 
erm is not divergent as « s -> 0, but it does vary 
with a s . We find that the linear T-dependence 
at small T is removed by —Aw 2 — » —Aw 2 + b\T, 



with b\ = 0.425 in dimensionless units. Next, 
we fixed the finite counterterm oc T 2 . We note 
that in the absence of any finite counterterm pro- 
portional to T 2 the second order phase transition 
(for H = 0) occurs at T c ~ 130 MeV, which is 
somewhat lower than that expected in QCD |||(| . 
We have pushed T c up to T c ~ 150 MeV by 
introducing -Aw 2 -> -\v 2 + b x T + b 2 T 2 , with 
&2 = —0.066 in dimensionless units. These values 
of bi and 62 were obtained with a lattice spacing 
a s = 0.4 on a 64 3 lattice, as shown by the filled 
symbols in Fig. [j]. On a 26 3 lattice with a s = 1 
(with the same physical volume) we find that with 
62 = 0.084 and b\ unchanged from above, the or- 
der parameter as a function of T is the same as 
that for a s = 0.4 within error bars (open sym- 
bols in Fig. 0). We have also verified that the 
a s -dependence of 62 vanishes for small A and is 
non-divergent in the a s — ► limit. 

Next we measured the T dependence of the 
masses and decay rates of the pion and a fields. 
These quantities follow from the dynamical re- 
sponse to small amplitude perturbations around 
equilibrium. Then the response takes the form 



b a (t — 0)) exp[— t/r] cos(wt), 



(5) 



with ui — y m 2 (T) — t~ 2 , for the particular case 
of a spatially homogeneous perturbation. Results 
are shown in Fig. |2| together with the chiral per- 
turbation theory prediction for the pion mass pf 



m 2 (T)=m 2 (0) 1 + 



T 2 



T « U . (6) 



The a becomes lighter in medium as a result 
of symmetry restoration at high T. m a has a 
minimum at T cross ~ 180 MeV, which defines the 
crossover temperature. At Tcross, fn-K — 220 MeV 
and m a ~ 280 MeV. The two masses are equal 
within error bars for T > 220 MeV. 

The timescale r(T) for the decay of a region 
of disoriented chiral condensate in contact with 
an equilibrated gas of hadrons has been com- 
puted previously by Steele and Koch using a 
hadron gas model and perturbatively (at 2-loops) 
by Rischke |J. We seek to choose rj(T) so that the 
t(T) we measure reproduces their calculations. 

A difficulty is that the classical thermal field 
theory already leads to long wave length decay. 
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Figure 2. Thermal masses of the a and tt fields vs. 
T. Error bars denote statistical and best fit uncer- 
tainties. The line shows the prediction from chi- 
ral perturbation theory, Eq. (|6|). The minimum of 
m a is reached at T C ross — 180 MeV. Calculations 
were done with a s = 0.4 on a 64 3 lattice. 

The associated time scale is much shorter (larger 
dissipation) than the target quantum r at low T, 
because of the absence of Bose-Einstein suppres- 
sion of short wave lengths in the classical case. 

Because 77 > in the Langevin equation (1) we 
find that our model can only be used to describe 
the long wavelength dynamics in the presence of 
a heat bath with T > 145 MeV, corresponding 
to about 80% of T cross . Above this temperature 
we can adjust our input bare rj to reproduce the 
DCC decay time computed by Steele and Koch, 
see Fig. |3|. However, T ~ 145 MeV is large enough 
that the assumptions in the calculation of r as in 
Ref. Q may be starting to break down, which 
means that our estimate of the limit of validity 
of our analysis may have a little play in it. The 
same considerations apply to the perturbative 2- 
loop calculation, with the additional caveat that 
at A — 20 it may be unreliable. 

In spite of these difficulties we are currently 
studying the non-equilibrium dynamics of the sys- 
tem induced by the combined effect of the cool- 
ing of the thermal bath and by volume expansion. 
The latter becomes the principal means of cooling 
at low T, which may help mitigate the problems 
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Figure 3. Circles show the DCC decay time r c i(T) 
for the classical a-model at 77 = 0. Squares show 
t(T) from the Langevin equation with rj{T) taken 
from Rischke's perturbative calculation. Dia- 
monds show t(T) from the Langevin equation 
with r)(T) chosen to reproduce t(T) from the cal- 
culation of Steele and Koch (solid line) . 

associated with the choice of r\. 
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